NONLINEAR INVERSION FROM PARTIAL EIT DATA: 
COMPUTATIONAL EXPERIMENTS 



S. J. HAMILTON AND S. SILTANEN 



Abstract. Electrical impedance tomography (EIT) is a non-invasive imaging method in which 
an unknown physical body is probed with electric currents applied on the boundary, and the 
internal conductivity distribution is recovered from the measured boundary voltage data. The 
reconstruction task is a nonlinear and ill-posed inverse problem, whose solution calls for special 
regularized algorithms, such as D-bar methods which are based on complex geometrical optics 
solutions (CGOs). In many applications of EIT, such as monitoring the heart and lungs of 
unconscious intensive care patients or locating the focus of an epileptic seizure, data acquisition 
on the entire boundary of the body is impractical, restricting the boundary area available for 
EIT measurements. An extension of the D-bar method to the case when data is collected only 
on a subset of the boundary is studied by computational simulation. The approach is based on 
solving a boundary integral equation for the traces of the CGOs using localized basis functions 
(Haar wavelets). The numerical evidence suggests that the D-bar method can be successfully 
applied to partial-boundary data in dimension two and that the traces of the partial data 
CGOs approximate the full data CGO solutions on the available portion of the boundary, for 
the necessary small k frequencies. 

1. Introduction 

1.1. EIT and the inverse conductivity problem. Electrical impedance tomography (EIT) is 
a non-invasive imaging method where an unknown physical body is probed with electric currents, 
and the internal distribution of conductivity is recovered from the measurement data. The 
reconstruction task is a nonlinear and ill-posed inverse problem, whose solution calls for special 
regularized algorithms, such as the D-bar method |31| . Applications of EIT include monitoring 
the heart and lungs of unconscious intensive care patients, industrial process monitoring and 
underground prospecting. 

Practical considerations typically restrict the boundary area available for EIT measurements: 
for example, it is not sensible to cover a patient completely with electrodes when imaging the 
heart. In this paper we study a possible extension of the D-bar method to the case when data 
is collected only on a subset of the boundary. See Figure [TJ The mathematical model of EIT is 
the inverse conductivity problem introduced by Calderon [ID]. Let f2 C M. n be a bounded and 
simply connected set with a smooth boundary d£l. Let a : 0, — > M be an essentially bounded 
measurable function satisfying cr(x) > c > for almost every x G VL. Let u £ H 1 ^) be the 
unique solution to 

(1) V • aVu = in f2, 

(2) u\ m = ^^(Sfi). 

The inverse conductivity problem is to recover the conductivity a from the Dirichlet-to-Neumann 
(D-N) map defined by 

A a : (p H> a— 

ov an 

Here v = (yi, v<i) = v\ + iv2 is the unit outward facing normal vector to the boundary. Here <ft is 
a voltage distribution we apply on the boundary, and A a (j) is the resulting current flux through 
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Figure 1. Left: original conductivity. Middle: reconstruction from full- 
boundary data using the D-bar method. Right: reconstruction from partial- 
boundary data using the proposed method. The subset T C d£l where the mea- 
surements are available is denoted by a black line, which corresponds to 25% of 
the entire boundary 

the boundary. Therefore, A a can be seen as an ideal infinite-precision model of practical voltage- 
to-current measurements. 

Calderon asked two questions in his seminal article |10j : 

(i) Is cr uniquely determined by A CT ? 

(ii) If the answer to (i) is yes, how can one calculate a from A CT ? 

In practical EIT imaging only a finite-range and noisy approximate operator is available. In 
general, A^ is not the D-N map of any conductivity. We usually only know that ||A^ — A,j||y < 5. 
Here Y is an appropriate data space and 5 > can be determined from the properties of the 
measurement device. This leads us to a third question: 

(iii) Given A^ and 5, how can one design a continuous map from Y to L°°(f2) whose output 
is a useful approximation to cr? 

As the inverse conductivity problem is ill-posed, the forward map *A : o i — \ A a does not have a 
continuous inverse. Therefore, question (iii) needs to be answered by constructing a regulariza- 
tion strategy [J3j. More precisely, a family of continuous mappings lZ a : Y —> L°°(f2) must be 
defined, parameterized by < a < oo, such that 

(3) lim \\K a (A a ) - a\\ La0 (Q) = 0, 

for each fixed a. Note that ([3]) is closely related to question (ii) above. Furthermore, one needs 
to specify a choice a = a{5) of regularization parameter as a function of the noise level so that 
a(S) — > as 5 — > 0. Finally, the reconstruction error \\lZ a ^(A^) — cr\\ L oo^ must vanish in the 
zero noise limit: for any fixed a we must have 

(4) sup {\\K a{s) (A 5 a ) - <7|| L oo ( n) : \\Al - A a \\y < 5} -> as 5 -> 0. 
AieY 

For more details, see Figure [2] and [3T| [37]. 

1.2. D-bar methods for full-boundary data. From the practical viewpoint, the solution of 
the inverse conductivity problem is a computational algorithm that implements a regularization 
strategy 1Z a satisfying ([3]) and Achieving such a goal is typically a large project involving 
several milestones, often of one of the following two types: 

(a) A theoretical breakthrough that outlines a computational approach 

(b) Successful computational experiments that inspire further theoretical study 

Let us review the history of a specific two-dimensional D-bar method for EIT in light of (a) and 
(b). 
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Figure 2. Schematic illustration of nonlinear regularization of the EIT problem. 
Here the forward map is defined as A(a) = A CT with domain of definition denoted 
by T>(A) C L°°(fi). The conductivity a is approximately recovered as 7Z a (A s a ). 

1996(a): Nachman showed uniqueness (i) and introduced a infinite-precision reconstruction 
method (ii) for twice differentiable conductivities in [40J. The proof used a nonlinear Fourier 
transform based on so-called complex geometrical optics (CGO) solutions, first defined by Fad- 
deev in 1966 [15J and later rediscovered in 1987 by Sylvester and Uhlmann in the context of 3D 
EIT |43| . Thus, [ID] represents a breakthrough in the form of (a) since it is the basis of the first 
numerical D-bar method |41[ l4"2l [36] . 

2004(b): Isaacson et al. demonstrated in |24[ [23] that the D-bar method performs well 
on practical data measured from laboratory phantoms and from human subjects. The manda- 
tory regularization step was provided by low-pass filtering in the nonlinear frequency domain. 
The need for such filtering is evident from the structure of the experimental nonlinear Fourier 
transforms: they blow up outside a disc centered at the origin. 

2009(a): The numerical evidence from practical imaging experiments [24\ [23] inspired a 
rigorous regularization proof of convergence in the form of ([!]), see |31] , This gave an answer 
to (iii) and outlined a method for choosing the regularization parameter as the inverse of the 
nonlinear cutoff frequency. We outline the reconstruction method in Section [2] below. 

There is an analogous history for other uniqueness proofs and related algorithms in two- 
dimensions. We review them briefly below without specifying explicitly progress steps of types 
(a) and (b). 

Brown and Uhlmann were able to prove uniqueness for real-valued conductivities assuming 
only one derivative in [9]. This result was complemented by constructive steps and numerical 
implementation by Knudsen and Tamasan |33[ [26] 127] : see also [32]. Francini [16J extended 
the uniqueness proof to complex conductivities whose real and imaginary parts are twice dif- 
ferentiable, and her approach was subsequently implemented in [171 [THJ, [19]. We outline this 
reconstruction method in Section [3] below. Both methods involve transforming |5| to a first 
order system of d z and d z equations. 

Astala and Paivaxinta answered Calderon's questions (i) and (ii) in their original smoothness 
category a 6 L°°(f2), see (4j[3]. This approach has been implemented numerically as well [2J[I]. 

Despite the above developments, some questions still remain open: 

• Is it possible to give a regularization analysis (iii) for less smooth conductivities than 
twice differentiable? There is numerical evidence of type (b) available since all of the 
above EIT methods produce noise-robust images when applied to data arising from 
discontinuous conductivities and regularized by nonlinear low-pass filtering [301 l29"j Wl\ [T] . 

• Can the D-bar methodology be used in the case of partial-boundary data? We discuss 
this in Section 11.31 below in the two-dimensional case. 
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1.3. Extension to partial-boundary data. It is of high practical importance to be able to 
compute EIT reconstructions from data measured only on a part of the boundary. One possibility 
for designing such algorithms would be to take one of the recent theoretical breakthroughs, such 
as |28 } 1251 [39[ I22j. and implement it in the spirit of (a) above. However, we do not discuss such 
approaches in this paper. We proceed along (b) and produce novel numerical evidence suggesting 
that it may be possible to use the classical D-bar approach for partial data reconstructions. It 
is our hope that these computational results inspire further theoretical advances. 

Our starting point is the assumption that only a proper subset T C d£l is available for 
measurements. We consider voltage-to-current data represented ideally by the restricted D-N 
map A a , defined as follows. Let <f> 6 H 1 / 2 (dCl) satisfy supp(0) C T and let u € be the 

unique solution of the conductivity equation 

(5) V • crVu = in fl, 

(6) 

Our partial D-N map is then defined by 

f »\ 7 7 du 

7 A a : ^ a— 

ov 

The practical data is a finite-range and noisy approximate operator satisfying || A^.— A CT ||y < 5. 

Let us briefly explain our approach in the context of the regularized D-bar method [31] based 
on Nachman's uniqueness proof |40| . In the full-boundary data case it begins by solving this 
Fredholm integral equation of the Second Kind for the (approximate) traces of the CGO solutions 
on dQ,: 

(8) 1>(z,k)\an = e ikg - I G k (z-0(^-MM(,k)dS(t), 

Jan 

where Gk is the Faddeev Green's function [15], here defined in the sense of tempered distribu- 
tions, 

(9) G k (z) :=e ikz g k (z), g k (z) := — 1 -= [ 6 dg. 

In the case of partial-boundary data, we solve the following equation for the unknown functions 
u(-,k) : T -)■ C: 



(10) u(z, k) = e ikz \ r - J G k {z - Q{1 5 C - A 1 )lo((, k) dS(Q, 
Now the hypothesis is that 

(11) ^(-,A;)|r ~ zercdtt. 



z g r c on . 



If (jlip holds, it opens up a variety of extensions of D-bar methods to partial-boundary data 
applications. 

1.4. Focus of this paper. How does one solve (jlOp numerically? Computational solution 
methods for boundary integral equations (BIEs) of type ©, corresponding to the continuum 
model, have most often been based on representing the unknown CGO solutions in terms of 
(generalized) trigonometric bases, where the basis functions are essentially supported on the 
entire boundary [311 [H [37]. This approach is not directly applicable to partial data problems. 
In this work we present new numerical experiments involving the solution of the above-mentioned 
BIEs using localized basis functions supported only on a subset of the boundary, in this case the 
Haar wavelets which are naturally applicable to the partial (as well as full) boundary continuum 
model. See Figure [3l 
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Trigonometric basis functions and 
corresponding electrode inputs: 
only for full-boundary data 
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Haar wavelet functions and 
corresponding electrode inputs: 
for both full and partial data 
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Figure 3. Illustration of various basis functions used for solving boundary in- 
tegral equations. Also shown are voltage patterns applied using 32 electrodes, 
approximating the basis functions. The right-hand side functions and patterns 
are localized and therefore more suitable for working with partial-boundary data. 



Let us stress that at present there is no proof available for the solvability of equation (jlOp . 
However, we did not encounter any problems when numerically solving (|10p . suggesting that it 
may be possible to prove unique solvability under appropriate assumptions. 

We demonstrate that it is possible to recover the traces of the CGO solutions approximately 
on the part of the boundary available for measurements. In other words, the approximation 
in (jlip is quite good in the C 2 and L°° conductivity examples we consider. In addition, we 
show below that these partial traces lead to interesting and useful reconstructions of practically 
relevant L°° conductivities. 

Our new results may be useful in extending three-dimensional D-bar reconstructions, such as 
[TTJ [HI [13], to partial-boundary data. 

We mention that numerical reconstructions using restricted information about the conductiv- 
ity have been published in cases of partial-boundary data, see [351 [ZD [2Ql S3] . The present work 
differs from those in that we aim to recover the full unknown conductivity function instead of 
inclusions in a known background. Also, there is an alternative methodology for partial-data 
EIT based on resistor networks [3i\ [7[ E]; our work again represents a very different approach. 
Finally, we mention that there is a large body of work on iterative solution methods for EIT in 
cases of partial data; those studies are fundamentally different from our direct (non-iterative) 
approach. 

2. Method 1: The method based on the Schrodinger equation 

As mentioned above, the proof by |40j transforms the conductivity equation ([5j) to the 
Schrodinger equation 

(—A + q(z)) ip(z, k) = 0, zeR 2 and A; G C, 
via the change of variables ifi(z, k) = yfou{z, k), where 



(*) 
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denotes the Schrodinger potential. Existence and uniqueness are then studied for the well known 
Schrodinger equation with CGO solutions tp that are asymptotic to e lkz and kz = (k± + ik 2 ){x + 

The alternative Lippmann-Schwinger formulation 

(12) n(z, k) = 1 - g k * (<?/u) , zeR 2 and k £ C, 

of the Schrodinger equation uses the related CGO solutions /i(z, k) = e~ tkz tp(z, k) where gk is 
related to the Faddeev Green's function and defined in ([9]). 

The reconstruction method of Nachman [JO] from infinite-precision data consists of the fol- 
lowing two steps: 

A CT — > t(k) — ->• a. 

Step 1: From boundary measurements A CT to the scattering transform t. For each fixed 
k £ C, solve in if 1//2 (<9f2) the integral equation 

(13) i>(z,k)\an = e ikz - I G k (z-0(K-AxMC,k)dS(C), 

J an 

for the CGO solutions V where the D-N map of the homogeneous conductivity 1 is 
denoted by Ai< Then, substitute ip into the formula for the nonlinear scattering transform 
t:C^C: 

Akz i 



(14) t(fc) = / e lk \K - AMz, k) dS{z), 

J an 

where dS denotes arclength measure on dCl. 
Step 2: From the scattering transform t to the conductivity a. Denote e(z,k) := 
exp(i(kz + kz)). For each fixed z£fi, solve the integral equation 

d5) ,( Z , k) = i + ^ I y)W^)dk' ldk ' 2 , 

then the conductivity is recovered by cr(z) = /i(z,0) 2 . 

The integral equation (|15p was obtained from a corresponding partial differential equation, a 
so-called D-bar equation, which involves the derivative with respect to the complex variable k. 
This is where the D-bar method gets its name. 

3. Method 2: The method based on a d z and d z system 

As the existence/uniqueness result by Francini [16] holds for complex admittivities 7 = a+iue, 
where cj is the frequency of the applied current and e denotes the electrical permittivity, and is 
an extension of that by Brown and Uhlmann [9] we will formulate the problem in the complex 
case. Let ui(z, k) and u 2 (z, k) be two CGO solutions to © with large k or z asymptotic behavior 
^rr and 

ik —ik 



£A and respectively. Introduce a matrix ^(z, k) of CGO solutions related to u\ and u 2 by 



*H ^ _ ^V2 ( 9 zUl \ / * u \ 1/2 / d z U 2 

*21 J 7 V d * U l J ' V ^22 J 7 V 

The transformed system is then 

(16) = QV, 

where D is a matrix of d z and d z partial derivatives and Q represents a matrix potential 
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Existence and uniqueness of solutions are then studied for (|16p instead of ©. As it is more 
practical to work with CGOs with finite asymptotic behavior, we often make use of the related 

matrix of CGO solutions M(z, k) ~ ( J ^ \ defined by 

(18) M(z k) - 9(z k) ° ^1 - ( e ~ iZk ^^ k ) e l ~ zk *Mz, k)\ 

(18) M[z,k)-V[z,k)^ Q e izk)-y e -i z kis, 2i{z k) eizk^fak)) ■ 

Similarly to Method 1, the full data direct reconstruction algorithm [17J also involves solving 
Fredholm integral equations of the Second Kind for CGO solutions using D-N data, evaluating a 
nonlinear scattering transform S(k), solving a equation, and using the recovered CGO solu- 
tions at k = to reconstruct the conductivity and permittivity. The method can be summarized 
in the following steps: 

A 7 S(k) A M{z,0) A 7. 



Step 1: From boundary measurements A 7 to the scattering transform S(k). 

For fixed k £ C\0, solve Fredholm integral equations of the Second Kind on d£l for the 
traces of the CGO solutions u\{z,k) and U2(z,k): 



(19) 
(20) 



Akz 



ui{z,k)\ d n = — 



an 



mi 



G k (z-0(lh-M)ui{C,k)dS{$ 



-ikz 



\an 



-ik 



(21) 
(22) 

(23) 
(24) 

Step 2: 

(25) 



G fe (-z + C)(A 7 -A 1 )« 2 (C,fc)d ( S(C). 

an J an 

Use the traces of u\ and 112 to compute the off diagonal entries of the CGO solutions 
$>(z, k) for z £ d£l from the BIEs 

gjfc(z-C) 



;>Q 4vr(z - C) 



[A 7 -A x ] u 2 (C,k) dS(() 



mi 



4ir(z - C) 



[A^-AiluiCC,*;) dS(Q, 



and compute the off-diagonal entries of the scattering matrix S(k) 



Si2(k) = — I e 
27r Ian 



S 2 i(k) 



-ikz 



ikz 



$i 2 0,A:)(i/i +iv 2 )dS(z) 
i $>2i(z, k)(vi - iu 2 )dS(z). 



27T J an 

Interpolate the scattering data S(k) to include k = 0. 

From the scattering transform S(k) to CGO solutions M(z, 0). 

Solve the equation (|25p for the matrix M(z, k) 



d k M(z,k) = M(z,k) 



e(z, k) 






e(z, —k) 



S(k). 



Step 3: From CGO solutions M(z, 0) to Admittivity 7 

Reconstruct the matrix potential Q from 

(26) Q12(Z)= M.(s,0) ' Q21[Z)= M + (z,0) ■ 
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where 

(27) M + (z,k) = M u (z,k) + e~^ kz+ ~ k ^M 12 (z,k) 

(28) M„{z,k) = M 22 (z,k) + e^ kz+k ^M 2 i(z,k), 
and use either Q± 2 or Q 2 ± to recover 7 

(29) 7(z) = exp <^ / — — T d/x(C) f = exp <^ / — dp(Q 

{ k J n z-C J I fin 2-C 

where the integration takes place over Q rather than all of C due to the compact support 
of the matrix potential Q. 

4. Computation of Partial Boundary Data CGO Solutions 

In this work we use localized basis functions in place of global basis functions. As mentioned 
above, the most commonly used global basis functions for the continuum model are the ex- 
ponential trigonometric basis functions e m6 '. When electrode models (such as the gap, shunt, 
or complete electrode model) are used, a trigonometric basis of sines and cosines is often used 
instead |12l 124] . As mentioned above, the common thread of these global basis functions is that 
their support is essentially the entire boundary of the domain. By contrast, localized basis func- 
tions are supported on a subset of the boundary. Examples of localized basis patterns include 
the skip patterns and adjacent patterns (see e.g., [IHl EU [19]) as well as the the Haar wavelets. 
In this work we use the Haar wavelets as they are localized basis functions that can be naturally 
used in both the continuum and electrode model cases. 

As the boundary integral equations in Methods 1 and 2 are very similar, we will describe, 
without loss of generality, the computation in detail for Method 1. In order to solve the boundary 
integral equation (fi~3"j) 

1>(z, k)\on = e lkz - [ G k (z - C) {K - Ai) V(C, k) dS(C), 
Jan 

for the traces of the CGO solutions ip(z, k), we will need the Dirichlet-to-Neumann (D-N) map, 
and thus we must first discuss the applied voltage patterns, in this case, the Haar wavelets. 

4.1. Description of Haar Wavelets. Let T denote a subset of the boundary d£l and let 
|r| = C denote the length of the subset T. The first wavelet is the constant wavelet which we 
will denote (fri, and is defined as: 

<fo.(z) = hi zircon 



(30) 



h 



1 - VI 



If z G dQ\T, the wavelet (f>x(z) is set to zero, as will be the case for the subsequent Haar 
wavelets. 

For ease of notation, let d(z) be the distance, along the subset T of the boundary, a point 
z is from the beginning point zq on T (corresponding to the smallest 9 value in the traditional 
counter-clockwise orientation) and zc the ending point (corresponding to the largest 9 value). 
Thus, d = at zq and d = C at Z£. 

The second wavelet is the so-called mother wavelet which we will denote 4> 2 and is defined as: 



(31) <h{z) 



hi, {z e r|o < d(z) < §} 

-hi, {z £ T\£ < d(z) < £} 



As the basis functions need to be orthonormal, we require ((j> m ,(f) n ) = o~m,m which the above 
wavelets satisfy by construction. 
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The third and fourth Haar wavelets 03 and 04 are copies of the mother wavelet 02, squished 
into 1/2 the length of the support of 02 as follows: 



»(*) 



(32) Mz) _ jh 2 , {z€r\$<d{z)<%} 




{z £ r|0 < d(z) < f } 

{ z er\±<d(z)<j} 



{z€T\f <d(z)<£}, 



h 2 = ^. 

Notice that these new wavelets satisfy {(f> m , <f>n) = 5m,n for m,n = 1, . . . , 4. 

The fifth through eighth Haar wavelets 05, 06 , 07, and 08 further subdivide the domain this 
time squishing into a region of support 1/2 that of 03 (for 05 and 06), and 04 (for 07 and 0s) 
as follows: 



(h 3 , {z G r|0 < d(z) < 
\-h 3 , {z£F\£ <d(z) < f }, 

{ z£ T\f<d(z)<f} 

{zev\% <d{z)<% }, 




(33) ffcs, {zer|f <d(*)<f£} 



07 (z 



-h 3 , { zG r|f <d(z) < f }, 
f/i 3) {^er|f <d(z)<f } 

{-h 3 , {zeT\ 7 ^ <d(z)<£}. 



Notice again that these Haar wavelets satisfy the orthonormality condition and from these 
formulas it becomes obvious how to extend the pattern to form the subsequent wavelets. An 
exact formula for the j'-th height function hj 



2J- 1 

corresponding to Haar wavelets with support width Wj 

c i>2. 



J 2^- ] 

4.2. Formation of the D-N map Using Haar Wavelets. As the main goal is to use only 
partial boundary data, thus applying and measuring data only on subset of the domain, it is 
more natural to apply voltages (rather than currents) and form the Dirichlet-to-Neumann (D-N) 
map directly. Although in practice currents are frequently applied, and thus the Neumann-to- 
Dirichlet (N-D) map is formed first (which is done to dampen noise), that approach requires 
the inversion of the N-D map, which for partial data poses new questions. As a preliminary 
approach, we proceed with Dirichlet data. 
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The conductivity equation ([5]) can be solved using the Finite Element Method. For each Haar 
wavelet, the Dirichlet boundary value problem is solved and the resulting solution u in Q is used 
to determine the current flux (Neumann data) at the boundary. This allows the determination 
of the Neumann data corresponding to the prescribed Dirichlet data and formation the D-N 
map: 

Kf = a 



d, 



v 



Note that for the cases in this document, the conductivity on and near the boundary is 1. 

The discrete matrix approximation to the D-N map is formed using the following formula for 
the (m, n)-th entry 

(34) A^(m,n) := (A CT c/> m ,0 n ) = (aV(/) m • v,<f) n ) = (V(p m • v,(j) n ), 

where 4>j are the Haar wavelets described in Section 14.11 that now serve as the Dirichlet data. 
As v denotes the outward facing unit normal and f2 is the unit disc, at the boundary point 
z = e %e = cos(#) -Msin(#), we have v = (cos(8) , sm(9)} . 

4.3. Solution of the Full Data BIE. After forming the D-N map using the Haar wavelets, 
localized basis functions, we proceed to solving the boundary integral equation for the traces of 
CGO solutions t/j(z, k). This involves the solution of a Fredholm integral equation of the Second 
Kind. Following the approach of [TH [T71 [13 LTB] , expand the exponential e lkz and CGOs ip(z, k) 
in the Haar wavelet patterns {^}^ =1 - Let be an evaluation point and J denote the number 
of linearly independent Haar wavelet functions. Then, 

(35) 

and 

J 

(36) e^wJ^-CAO^, 

where $ denotes the normalized Haar wavelets such that their £ 2 norm is 1, i.e. they are related 
via 

Hi 



Let h(k) denote the column vector h(k) = [bi(k), . . . , bj(k)] T , and define c(k) analogously where 
T denotes the standard matrix non-conjugate transpose. 

Let Ei> denote the l'-th subdivision of the boundary <9$7 (£' = 1, . . . , L) centered at the center 
of the ^'-th boundary element z^ with length 2ir/L. Splitting the integral over d£l into a sum 
of integrals over the subsections Eg 

1>(zt, k) « e lkz * -J2 [ G k (z e - C) SKMCt' , k) dS(Q 
(.1=1 E i' 

= e ikz t _ I G k { Zl - C) dS(() [5A<r1p(Qr,k)] , 

£t_i J Erf 



where for ease of notation 



SA a = K - Ai. 
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Using the expansions for ip(z£,k) and e tkze , ([35]) and ([36]) respectively, we have 



3=1 3=1 



= 1 



3=1 



E<y(*o*<-E / G *(*<-o ^(c)E^(^(o 

7=1 e=i jE i' 7=1 



where /j (^') denotes the action of the discretized <5A^ matrix on the j-th normalized Haar 
wavelet basis function evaluated at Q>. Define the matrix approximation to the Faddeev Green's 
function as 



(37) 



G k (£,£') 



G k (zt,Ce) i±e 
£ = £' 



removing the singularity at Gfc(O). Then 
J J 



(38) 

Following [12] 
(39) 



Itt 



3=1 



3=1 



3=1 



£'=1 



/p(60 « (MA**) 



i.e., the (£', j) entry in the matrix resulting from multiplication of the matrix of normalized basis 
functions $ and the discretized difference in D-N maps SA^f . Using the properties of matrix 
multiplication, equation (j38[) can be rewritten as 



2vr 



Em*o*j = E c ;(*o^ - -rE w G fc ^ A * f 



3=1 



3=1 



L 



3=1 



or equivalently, 



2vr 



$b = $c - ^G fc $5A^b, 



a matrix equation for the unknown coefficients b which are needed in the normalized Haar 
wavelet basis expansion of ip(z, k). 

Using the orthonormality of the normalized Haar wavelet basis functions in the matrix <I>, we 
multiply both sides of the equation by <& T , and then solve 



(40) 

where 

(41) 



(I + A)b = c, 



2ir 



To reiterate, for each desired value of k 6 C, expand e lkz for z £ dQ. in the normalized Haar 
wavelets to define the vector of coefficients c, and solve the system (I40p using GMRES for the 
unknown coefficients b. These coefficients are then used to reconstruct ift(z, k) for the specified 
value of k via (1351). 
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4.4. Solution of the Partial Data BIE. We now proceed to the problem of interest, namely, 
the solution of the boundary integral equation (|13p when only part of the boundary is accessible 
for data acquisition. The solution method is nearly identical to the full data Haar wavelet case 
presented above. 

Now r is a proper subset of the boundary d£l and the Haar wavelets and D-N map are formed 
as above. The D-N map now corresponds to data taken only on the proper subset T since the 
applied voltage is off T. 

Let z denote the boundary values z restricted to T and if) the corresponding partial data CGO 
solutions. We then expand e tkz and tp as before and solve the resulting system for each desired 
value of k G C: 

(/ + i)b = c , 

where 

(42) A = j$ T G k $5AM, 

and C = |r|. 



5. Numerical Reconstruction of Conductivities from Partial Data Using 

Method 2 

As stated above, the boundary integral equations (fl~9j) and (f20|) are nearly identical to the 
BIE (fl~3j) described above for Method 1. Thus, their traces on the subset T of the boundary 
can be recovered by solving analogous formulas. The coefficients b 1 for u\ are determined by 
solving 

(I + A)b 1 = c 1 , 

where c 1 are the coefficients in the Haar expansion of ^nr- Similarly, the coefficients b 2 for U2 
are determined by solving 

(I + l 2 )b 2 = c 2 , 

where c 2 are the coefficients in the Haar expansion of and A2 now contains the matrix 

approximation of Gf.(—z + C) instead of Gk{z — C)- 

A natural question is whether these CGO solutions, which match very well on V with their 
full data counterparts, can be used to produce informative reconstructions of the conductivity 
(and/or permittivity) near the region of the accessible boundary. Our aim was to understand the 
extent of the impact of the partial data CGO solutions on the remainder of a D-bar algorithm. 
Therefore, as an initial test, we left the remainder of the algorithm for Method 2 in tact which 
means computing the intermediate CGO solutions and ^21 using the partial D-N map and 
the partial data CGO solutions u\ and 112, computing the scattering transforms S12 and S21 
over r, and proceeding with Steps 2-3 as before. For the numerical details regarding how to 
implement Steps 2-3 of Method 2 see [18] . 

6. Computational Experiments 

We considered two test problems. Test 1 aims to determine how a partial data D-N map 
affects the values of the traces of the CGO solutions on the accessible portion of the boundary. 
Test 2 aims to determine the effect of the partial data CGO solutions on a D-bar algorithm. 
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6.1. Test 1: Partial Data Traces of CGO Solutions. For the first test problem we consid- 
ered the C 2 smooth conductivity given in Figure HI The conductivity equation was first solved 
using the Finite Element method for 256 Haar wavelets with essential support on the entire 
boundary, serving as 256 different Dirichlet boundary conditions. We considered the 3/4, 1/2, 
and 1/4 data problems with 192, 128, and 64 Haar wavelets respectively. Each of the partial 
data cases is centered around z = 1, i.e. 9 = 0. 

Using Method 1, we solved the full data matrix formulation (|40p of the boundary integral 
equation for the traces of the CGO solutions ip. The partial data traces of the CGO solutions were 
recovered by solving (|42p . In order to evaluate how well the reconstructed traces compare to the 
true traces we also solved the Lippmann-Schwinger equation (|12p using the twice-differentiable 
conductivity in Figure 01 

Figures [6] and [7] show the recovered traces of the CGO solutions if)(z,k) for k = 0.5, 4, 
and — M respectively, plotted against the true traces produced via the Lippmann-Schwinger 
computation. Preliminary results suggest that for small magnitude k, the partial data CGO 
solutions agree with the full data (and true) solutions on the accessible part of the domain. As 
the magnitude of the frequency parameter k increases, the partial data CGO solutions begin to 
drift slightly from the full data solutions. However, in the nonlinear CGO approaches typically 
used in EIT/EIS imaging, only low frequency CGO solutions are used and therefore these results 
are very promising. 




6.2. Test 2: D-Bar Reconstructions of Conductivities using Partial Data CGO Solu- 
tions. In Test 1, we saw that the traces of the CGO solutions computed using the partial D-N 
data corresponding to an accessible subset T of the boundary d£l match the true, as well as full 
data, traces very well on T for low frequencies k and C 2 smooth conductivities. Real life situa- 
tions often involve cases where the conductivities are not smooth but instead only bounded. In 
order to determine if the traces of the partial data CGO solutions still provide useful information 
when the smoothness is relaxed, we consider the L°° test conductivity shown in Figure El This 
test phantom could represent a saline filled tank containing an object of higher conductivity. 

As the smoothness condition for Method 2 is violated (as well as that for Method 1), we cannot 
compare the full or partial data traces of the CGO solutions to their "true" traces. Precisely, 
we cannot compute the potential Q (or q) and the associated d z — d z system (or Lippmann- 
Schwinger equation) for the "true" traces of the CGO solutions. Instead, we compare the partial 
D-N data traces of the CGO solutions to the corresponding full D-N data traces. Figures [9] and 
[TOl shows the reconstructed traces of the CGO solutions u\ and «2> respectively, for k = 3 + 2>i 
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Real parts of tp Imaginary parts of tp 




Figure 5. Traces of the CGO solutions ip resulting from various computations. 
Here k = 0.5. 

resulting from the conductivity distribution in Figure [8] plotted for full, 3/4, 1/2, and 1/4 D-N 
data with 256, 192, 128, and 64 Haar wavelets respectively. As the boundary integral equation 
for Method 1 is very similar to that of ui, we do not show it here. Clearly the partial D-N data 
traces of the CGO solutions appear to match the full D-N data traces of the CGO solutions in 
the accessible region T of the boundary. 

Next we used the partial D-N data traces of the CGO solutions u\ and U2 in the modified 
D-bar algorithm for Method 2, described above, to determine their effect on the algorithm and 
thus the reconstructed conductivity distribution. Figures [TT] and [12] show the reconstructed 
conductivity from full, 3/4, 1/2, and 1/4 D-N data using scattering data satisfying \k\ <= 3 
and \k\ <= 4, respectively. The range of the reconstructed values decreases with the size of the 
accessible region of the boundary T, however an object is clearly visible in all cases. As the 
magnitude of k increases the reconstructed values of the conductivity improve, however, as in 
the full data D-N case, increasing the scattering radius too much can introduce artifacts into 
the reconstruction. 

Note that the reconstructions are all plotted on their own scales. Viewing the reconstructions 
in this manner allows us to view the conductive object more clearly. From the figures we are 
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Figure 6. Traces of the CGO solutions ip resulting from various computations. 
Here k = 4. 



able to determine whether the object is more or less conductive than the background as well as 
its approximate location even from as little as 25% D-N data. 



7. Conclusion 

The first step in nonlinear EIT/EIS imaging uses the data to determine the traces of the 
CGO solutions at the boundary. This is done by solving a boundary integral equation which is 
a Fredholm equation of the second kind, e.g., (fl~3|h 

In this work we use simulated partial data and a wavelet-based integral equation solver to 
demonstrate that CGO solutions can be approximately recovered from partial data on the appro- 
priate part of the boundary. This result is clearly seen in Figures El [6] and [7] for a C 2 conductivity 
and in Figures l9l and flOl for an L°° conductivity. In addition, we have demonstrated that such 
partial data CGO solutions can be used in existing full data D-bar methods to provide useful 
and informative reconstructions, even in the case of discontinuous conductivities. 
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Figure 7. Traces of the CGO solutions tp resulting from various computations. 
Here k = —4i. 
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